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We discuss aspects of a covariant QCD modeling of meson physics by illustrating 
applications to several coupling constants and form factors. In particular, we cover 
the pTTn and 7r 77 interactions, the p contribution to the pion charge radius, and 
ttNN coupling. 



1 Introduction 

The Dyson-Schwinger equation (DSE) approacrE to non-perturbative QCD 
modeling provides semi-phenomenological gluon and quark two-point functions 
that have proved toJpe quite efficient in describing and predicting the physics 
of low-mass mesons 0. To outline the basis of such efforts, consider the fully- 
dressed and renormalized quark propagator defined by the Dyson-Schwinger 
equation (DSE) in Euclidean metric 

S- 1 (p)=Z 2 [i 7 .p+m (A)] + Z 1 ^ J* ^fD^(p-kyypS(k)r«(k,p), (1) 

where mo (A) is the bare mass parameter and A characterizes the regulariza- 
tion mass scale. Here the dressed gluon propagator D^(q) and the dressed 
quark-gluon vertex T 9 ^{k,p) are the renormalized quantities and they satisfy 
their own DSEs which require at least some higher order n-point functions. 
The ultraviolet behavior of D^ u {q) and r^(/c,p) can be constrained by pertur- 
bation theory. However physics at the hadron length scale depends crucially 
upon the behavior of S(p), D flu (q) and T 9 ^(k,p) in the infrared. Present QCD 
modeling begins with an explicit solution of the quark DSE via phenomenolog- 
ical IR forms for D flv (q) and T^(k,p). Although the DSE itself can probe the 
pseudoscalar component of the pion via the chiral condensate and dynamical 
chiral symmetry breaking, the Bethe-Salpeter equation (BSE) is required to 
address other mesons and to obtain the sub- leading pion components. 

The BSE for a bound state of a quark of flavor fi and an antiquark of 
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flavor /2 is 

r(p; P) = y ^4 #(p, g; P)S fl (g + £P)r(g; P)% (g - £P) , (2) 

where £ + £ = 1 describes momentum sharing. The kernel K operates in the 
direct product space of color, flavor and Dirac spin for the quark and antiquark 
and is the renormalized, amputated qq scattering kernel that is irreducible with 
respect to a pair of qq lines. At the present stage of QCD modeling, the BSE 
is employed in ladder approximation with bare vertices, that is 

K(p,q;P) = -g 2 D flv (p-q)^ ltl ® ^-j v . (3) 

The treatment of the quark DSE that is dynamically matched to this is the bare 
vertex or rainbow approximation for then the axial vector Ward-Takahashi 
identity is preserved and the Goldstone theorem is manifest. The various 
models of this type generally use the Ansatz 

where a e g(q 2 ) implements the pQCD running coupling in the UV and a phe- 
nomenological enhancement in the IR. The pion and kaon are well described 
in a recent work of this type.13 

2 QCD Modeling of Mesons and Interactions 

In Euclidean metric, the mass-shell condition for meson couplings requires the 
quark propagators in loop calculations be evaluated at complex quark momen- 
tum. To facilitate a broad survey of such applications, the present approach 
is to make use of a convenient analytic parameterization of confined solutions 
of the quark DSE. The broad features are taken from the solution to a simple 
DSE modclo that is extremely infrared dominant, produces a propagator with 
no mass-shell pole, and includes gluon-quark vertex dressing determined by 
the Ward identity. The resulting propagator is an entire function in the com- 
plex p 2 -plane describing absolutely confinedu dressed quarks in the presence of 
both explicit and dynamical breaking of chiral symmetry. The finer details of 
more realistic DSE solutions are accommodated by typically five parameters 
that are are used to restore a good description of pion and kaon observables: 
fn/K] m Tr/Kf\iffQ)] r ir] the 7r-7r scattering lengths; and the electromagnetic pion 
form factor. UU 
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The general form of the pion Bethe-Salpeter (BS) amplitude is 



Tl(k; P) = r J 75 lE^k; P)+y PF„(k; P) 



+ 1 -kk-PG„(k;P) + a fiV H n (k; P) 



(5) 



and the first three terms are significant in realistic model solutions @ and are 
necessary to satisfy the axial Ward-Takahashi identity, a The latter identity, to 
lowest order in P at the chiral limit, yields E„(k; P s= 0) = B(k 2 )/ where B is 
the scalar part of the quark dynamical self-energy. U It has been quite common 
to assume that only this term of is important for pion coupling. This 
has been questioned by recent studies of interactions, such as that shown in 
Fig. |l|, where we employ approximate it BS amplitudes such as those obtained 
from a rank-2 separable ansatza for the ladder/rainbow kernel of the DSE and 
BSE. They preserve Goldstone's theorem and should be adequate for infrared 
integrated quantities. Parameters are fit to m^/n and f^/K- The resulting 7r 
BS amplitude is 



r£(fc; Q) = klg{k 2 )\{ + i^f(k 2 )\ p 2 + i^e^xp^kxQMk 2 )^. (7) 




(G) 



The BS amplitudes are normalized in the canonical way. 




k+P/2-Q/2 



P-Q/2 



k-P/2 



Figure 1: Diagram for the pirn calculation. 
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3 The pirn, 7r°77 and jtttt Interactions 

The first term in a skeleton graph expansion of the pirir vertex can be 
expressed as 

a„(p, Q) = 2iv c tr s J ^s( ? ')r;s(,")r I s(g'")r, , (8) 

The momentum notation can easily be deduced from the Feynman rules asso- 
ciated with Fig. |, e.g. q' = k + P/2 + Q/2 and r£ represents T^(fc + P/2; Q) 
etc. With both pions on the mass-shell, P • Q = and P 2 =— to 2 . — Q 2 /4. 
In this case symmetries require the form A At (P, Q) =—P IJ ,Fp 7T7T (Q 2 ) and the 
coupling constant is <7 p7r 7r =Pp7r7r(Q 2 = — w 2 ). 

Previous investigations of the /97T7T coupling constant in terms of a covariant 
quark-gluon phenomenology for the intrinsic properties of p. and 7r employed 
only 7 M and 75 covariants for the respective BS amplitudes. LL30 It has since 
been demonstrated for a number of infra-red sensitive quantities such as to w 
and /jr, that the pseudovector terms Wr-the pion BS amplitude are responsible 
for corrections in the 20-30% range. mM The model p amplitude in Eq. (0) 
is admittedly crude, but the relative magnitude of the three surviving scalar 
amplitudes will hopefuly provide qualitative guidance. With the separable 

Table 1: g p7 nr calculation and contributions from meson covariants. 



Spir-ir = 6.28 [expt 6.05] 



7r Covariants 


p Covariants 


75 171% 
757 ■ Q -71% 


7m 94-5% 
75Cm 1 k Q 5.5% 
fc M 0.01% 



model BS amplitudes of Eqs. m) and (|7|), the prediction for g P inr, given in 
Table [l], compares favorably with the empirical value associated with the p — > 
tttt decay width. Truncation to the dominant p amplitude is found to only 
make a 5% error. However the sub-dominant pion component (pseudovector 
757 • Q) enters quadratically here and makes, a-major contribution (-71%). 

Studies of pion loops in the p — to sector EjtIj and in the pion charge form 
factorEj suggest that the qq extended structure of the pion significantly weakens 
such contributions compared to models or effective field theories built on point 
coupling. A similar issue arises in the role of the p in the space-like pion charge 
form factor. The dressed photon-quark vertex T u {q\ Q) can be separated (non- 
uniquely) into a p pole or resonant piece (which is transverse) and a background 
or non-resonant piece (which is both longitudinal and transverse). Thus in the 
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present approach the pion charge form factor takes the form 

where Hip(Q 2 ) is the fry polarization tensor and F^ IA (Q 2 ) is the generalized 
impulse approximation (GIA) result due to the non-resonant photon-quark 
coupling. It has been found to be phenomenologically successful in the space- 
like region and a persistent result is that 85 — 90% of the charge radius is 
naturally explained that way. El 

Is the p contribution small enough in the present QCD-modeling approach? 
WithF^(Q 2 ) = g^fp.AQ 2 ) andn£ 7 (Q 2 ) = -Q 2 f P1 (Q 2 ) / g v , which is con- 
sistent with electromagnetic gauge invariance, both form factors / depend on 
meson substructure dynamics and have been calculated. The p contribution 
to rv from Eq. (^) is then 

(rT' e ) 2 = rl - (r° IA ) 2 = 1.2 (0)^(0) ± , (10) 

m p 

where wc have used the empirics! result g /gy ~ 1.2 rather than univer- 
sal vector coupling. Our calculations include the extended nature of the 
mesons and produce / p7r7r (0) ~ 0.5 and / P7 (0) «s 0.65. This yields (rP ole ) 2 = 
0.16 fm 2 . In contrast, the empirical Vector Meson Dominance (VMD) picture 
has r^! IA — and f p7T7T = f pl = 1 so that r 2 ~ 6g„ 7r7r /(m 2 (7y) and produces 
~ 0.4 fm 2 . Adding the non-resonant impulse results (r^ IA ) 2 — 0.31 fm 2 gives 
a total of 0.47 fm 2 with our present approach. This is obviously an overesti- 
mate of the experimental value (0.44 fm 2 ) leaving no room for the pion loop 
contribution of the expectedE3 size. However, the main point is that a p con- 
tribution to the pion charge radius which is much smaller than that from the 
simple VMD assumption is consistent with the present status of DSE-based 
QCD modeling of the pion. 

The coupling constant for the tt° — > 77 decay is given by the axial anomaly 
and is a consequence of gauge invariance and chiral symmetry in quantum 
field theory. A modeling of nonperturbative QCD should preserve these fea- 
tures. This was the aim of a study Ej that applied the present approach and 
also investigated the form factor for the transition 7*71-° — * 7. The calculation 
employed a dressed quark loop similar to Fig. ^. The dressed photon-quark 
vertex was represented by the Ball-Chiu Ansatz which satisfies the relevant 
symmetries and obeys the Ward-Takahashi identity and is completely speci- 
fied by the amplitudes of the dressed quark propagator. Only the first term 
of Eq. (|J) for r w was retained and the chiral limit axial Ward identity result 
E^(k;P) -> B(k 2 )/f v was used. 
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It was verified that the axial anomaly result for the coupling constant was 
obtained (independent of the details of the quark propagator parameterization 
used) as a test of the numerical work. The remaining pion BS amplitudes 
from Eq. (||) evidently do not contribute to the chiral limit coupling constant 
and this has been recently demonstrated analytically. Ii.3 However such terms 
can become increasingly important at higher mass scales and, for example, are 
crucial for the asymptotic behavior of the pion charge form factor. liil For this 
reason, the large momentum behavior of the calculationlij of the axial anomaly 
transition form factor shown in Fig. ^| can be expected to receive significant 
corrections when the sub-leading pion BS amplitudes are included. This is a 
difficult task presently under study. 



4 irNN Coupling 

One can ask whether the dynamical content of the simple pion BS amplitude of 
Eq. (||) produces a irNN coupling constant consistent with the well-established 
empirical result g-nNN — 13.4. We make an estimate, using the valence quark 
states of a mean field chiral quark-meson model Ii3ii3 of the nucleon in which 
the chiral meson modes are generated as qq correlations. Thi^ approach has 
previously proved fruitful for the pNN and uiNN couplings. Ell In Euclidean 
metric, the ttNN vertex is 

K nNN (Q) = ±-(N\ J ^q(p+Q)TAp;Q)q( P -%)\N) (11) 
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where q is the quark field, Q is the tt momentum, \N) is the static mean field 
nucleon state and T T is the BS amplitude. The nucleon valence quark wave 
function renormalization constant Zn arises from the dynamical nature of the 
quark self-energy. Ej At the it mass-shell, rV is normalized in the canonical way 
such that it is the residue of the pseudoscalar qq propagator there. 

The standard form factor F^NNiQ 2 ) is identified from recasting the results 



from Eq. (11) into the form 



Knn{Q) = «Jv(f )[ i l5 T N F vNN (Q 2 ) ]« W (y) . (12) 
The nucleon mass shell condition does not allow distinction to be made be- 
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Figure 3: Pion-nucleon form factor normalized such that g^NN i s t ne value at the pion mass 
shell value Q 2 xs 0. The solid line includes both PS and PV components of the pion, while 
the dashed line includes the PS part only. 



tween pseudoscalar or pseudovector coupling, and only one form factor can be 
identified. Such is not the case for the constituent quarks and there will be 
distinct contributions to F-xnn from the PS and PV components of the it BS 
amplitude. As a mean field nucleon model cannot properly address the recoil 
issue, we resort to the common prescription in which Brcit frame kinematics 
is used to define a form factor a t lo w momentum transfer. 

In Fig. ||we show the resultlia for Ft T nn(Q 2 ) for space-like tt momentum. 
With both PS and PV components of the pion included we obtain g^NN ~ 11, 
while use of just the PS pion gives g^NN ~ 17. Rather than compare directly 
to the empirical value 13.4, the only conclusion we draw from this estimate 
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is that, without the PV component of the ir BS amplitude, g^NN would be 
overestimated by almost 50%. 

5 Summary 

Since the parameters in this approach have been previously fixed through the 
requirement that soft chiral quantities such as m^iKi fw/K an d charge radii 
r Tr/K be reproduced, the meson couplings discussed here have been produced 
without adjusting parameters. The results imply that this present approach 
to modeling QCD for low-energy hadron physics can capture the dominant 
infrared physics. The PV component of the pion is found to be important 
(at the level of about 25% and above) for a variety of physical quantities such 
as m T , /Vr, 3 pw7r , and g^NN ■ We expect that the large momentum behavior of 
form factors such as "/tttt and 7777 will require attention to both types of PV 
amplitude evident in Eq. (||). 
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